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ABSTRACT 

We present theoretical delay times and rates of thermonuclear explosions that 
are thought to produce Type la supernovae. including the double-detonation sub- 
Chandrasekhar mass model, using the population synthesis binary evolution code 
StarTrack. If detonations of sub-Chandrasekhar mass carbon-oxygen white dwarfs 
following a detonation in an accumulated layer of helium on the white dwarf's surface 
("double-detonation" models) arc able to produce thermonuclear explosions which are 
characteristically similar to those of SNe la, then these sub-Chandrasekhar mass ex- 
plosions may account for at least some substantial fraction of the observed SN la rate. 
Regardless of whether all double-detonations look like 'normal' SNe la, in any case the 
explosions are expected to be bright and thus potentially detectable. Additionally, we 
find that the delay time distribution of double-detonation sub-Chandrasekhar mass 
SNe la can be divided into two distinct formation channels: the 'prompt' helium-star 
channel with delay times <500 Myr (13% of all sub-Chandras), and the 'delayed' dou- 
ble white dwarf channel, with delay times >800 Myr spanning up to a Hubble time 
(87%). These findings coincide with recent observationally-derived delay time distri- 
butions which have revealed that a large number of SNe la are prompt with delay 
times <500 Myr, while a significant fraction also have delay times spanning ~ 1 Gyr 
to a Hubble time. 
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1 INTRODUCTION 

The exact nature of the stars that produce Type la su- 
pernovae (SNe la) - which are believed to be thermonu- 
clear explosions of carbon-oxygen (CO) white dwarfs (WDs) 
close to the Chandrasek har mass limit - remains unknown 
(e.g., iBranch et ail 1 19951 ). The most widely favoured SN la 
progen itor scenarios involve the double degen erate scenario 
(DDS; IWebbink|[l983 : llben fc Tutukovlll984r) and the sin- 
gle degenerate scenario (SDS: IWhelan fc Iberj|l973h . In the 
DDS, the merger of two CO WDs with a total mass ex- 
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ceeding the Chandrasekhar mass limit, Mch ~ 1-4 Mq, 
can lead to explosive carbon burning which causes a SN 
la explosion. In the SDS, a CO white dwarf accretes from 
a hydrogen-rich stellar companion via stable Roche-Lobe 
overflow (RLOF) and undergoes hydrogen burning on the 
surface, enabling the WD to accumulate mass toward Mch 
until carbon is ignited explosively in the centre of the WD 
leading to a SN la. However, in the stable RLOF configu- 
ration the companion does not have to be a hydrogen-rich 
main sequence (MS) or giant-like star, but can be a non- 
or semi-degenerate heliu m- burning star, o r a (degenerate) 
helium white dwarf (e.g., Iben et all 1 19871 : lYoon fc LangeJ 
120031 ; ISolheim fc Yungelsonll2005h . Like the SDS, the WD 
explodes once it approaches Mch (this helium-rich donor 
scenario will hereafter be referred to as HeRS). 

Recently, Ruiter et al. 2009 (Paper I) carried out a pop- 
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ulation synthesis study showing rates and delay times - time 
from birth of a progenitor system in a short burst of star 
formation to SN - for three formation channels of SNe la: 
DDS, SDS and HeRS. Additionally it is worth asking the fol- 
lowing question: to what degree do sub-Chandrasekhar mass 
WDs contribute to the population of explosive (SN la-like) 
events? We point out that in some scenarios, a layer of ac- 
cumulated helium on an accreting, sub-Chandrasekhar mass 
WD's surface may undergo sever al shell flashes (e.g., lead- 
ing up to a ".la" as discussed in iBildsten et al J 120071 ). and 
while these explosions can be bright, none of them result in a 
SN la. However, it is possible that the detonation in the he- 
lium layer causes the underlying CO WD to detonate, result- 
ing in a (final) explosion that does look like a SN la. Here, 
we extend our study of progenitors and focus our analysis 
and discussion on these thermonuclear "double-detonation" 
events involving helium-rich donors, hereafter referred to as 
the sub-Mch model, in which a CO WD accretes from a 
helium-rich companion filling its Roche-Lobe and explodes 
as a SN la before reaching the Mch limit. 

The sub-Mch model has thus far been regarded as an 
unlikely model for SNe la owing to the fact that most syn- 
thetic light curves and spectra of these objects from previous 
studies did not match those observed for SNe la. However, 
it has recently been argued that the model might be capable 
of producing a better match to observation, depending on 
details reg arding the manner in which the accre ted helium 
burns (e.g. iFink et~ai]|2010l : Ikromer et alj|20lj ). In either 
case, the explosion mechanism is expected to produce events 
that are bright and should be detectable. Thus, quantifica- 
tion of their predicted rates and delay times is an important 
step for testing our population synthesis models, and for 
determining what fraction of SNe la could conceivably be 
associated with this channel. 

Since these calculations are based on the work that was 
performed for Paper I, the reader is referred to that paper for 
a more detailed description of the DDS, SDS and HeRS sce- 
nario. The layout of this paper is as follows: In section 2 we 
summarise some background information on SN la progeni- 
tors from the literature. In section 3 we discuss the popula- 
tion synthesis modelling. In section 4 we present delay time 
distributions and rates as a function of stellar mass as well as 
distributions showing the exploding CO WD core mass. In 
section 5 we close with a discussion of these findings and pos- 
sible implications/predictions for SN progenitors and their 
host stellar populations. 



2 BACKGROUND 

2.1 Recent observations of SNe la delay times 

The idea that SN la pr ogenitors belong to at least 
two distinct populations (|Scannapieco fc Bildstenl 120051 : 



apiL _ . 

M annucci et al.ll2006l : |Pritchet et aI.ll2008T l has been gaining 
ground. A picture is emerging which supports populations of 
both quickly-evolving (prompt) progenitors with short delay 
times less than ~500 Myr, as well as more slowly-evolving 
progenitors with (sometimes rather) lo ng delay times span- 
ning up to a Hubble time (but see aIso lGreggidl201oT ). 

The delay time distribution (DTD) is a useful tool 
in determining the age of the progenitor stellar popula- 



tion, which places strong constraints on the different pro- 
posed progenitor scenarios. There are a growing number 
of observationally-derived DTDs present ed in the litera- 
ture from various groups (see section 1 of lMaoz fc Badenes 
|2010|. for an overview of these previous studies'). iTotani et al. 
(j2008f ) derived the DTD from a large population of old 
galaxies which only probed delay times >100 Myr, and 
they found that the DTD follows a relatively smooth 
power- law distribution (t -1 ) from ~0.1 to 8 Gyr (see also 

Pr obing younger stellar popula 



iHoriuchi fc Beacomll2010h. 

Maoz fc Badenesl (|2010h were able to determine that 



tions, 



a substantial fraction of SNe la are promptQ Among these 
prompt SNe la (35-330 Myr delay times in that study) the 
SN la rat«0is ~0.09-0.40 SNuM, compatible with the results 



of iLi et all (|2010al ) , whereas delayed SNe la in that study 
(330 Myr-14 Gyr delay times) had an overall smaller rate: 
<0.0024 SNuM. This study confirmed that roughly half of 
SNe la occur with delay times <330 Myr, thus g iving strong 
suppor t for a prompt component of the DTD. iMaoz et alj 
(|2010bl ) reconstructed the star formation histories for a sam- 
ple of LOSS SN host galaxies and found strong evidence for 
both a prompt la component with delay times < 420 Myr 
and a delayed comp onent with long delay times (> 2.4 Gyr). 
iBrandt" et all (|2010l ) used SN light -curves and spectra from 
host galaxies of 101 SNe la with 2 < 0.3 to construct the 
DTD, and arrived at a similar conclusion: that roughly half 
of SNe la occur with delay times <400 Myr, while the other 
half have long (> 2.4 Gyr) delay times. Further, they find 
that the short delay time events are more luminous with 
slowly-declining light-curves, and are associated with young 
stellar populations, whereas the SNe with long delay times 
are typically fast-declining, sub-luminous events. 

As one can see, the aforementioned recent studies in- 
dicate that SNe la are observed to occur over a range of 
long delay times (with less events as time goes on), with a 
substantial fraction also occurring at very early times. 



2.2 Two progenitor scenarios: DDS and SDS 



For 
tions 



s ome time, p o pulat io n synthesis c alcula- 



(llben fc Tutukovl Il984l : 



Yungelson fc Tutukovl 19971: lYungelson fc Livid 



Yungelson et all [l994; 

1998 



Nelemans et all 2001 : iRuiter et al. 20091 ) have predicted 



that the number of merging CO WDs with a total mass 
exceeding Mch (DDS) is sufficient to match, and thus 
possibly account for, the rate of SNe la (0.4 ± 0.2 per 
century for the Galaxy, ICappellaro et all I1999T ). At the 
same time, the theoretically-predicted SN la rate from the 
SDS channel is usual ly unable to explain the ob served rates 
of SNe la (see also iGilfanov fc Bogdanl |2010| ). There are 
very few SNe la that show any hint of hydrogen lines in 
their spectra; if the progenitor involved a hydrogen-rich 
companion, in particular a giant donor, one may expect 



1 A co mplimentary result was also determined by IMaoz et al.l 
feOlOal) , who found that in galaxy clusters the DTD is well-fit by 
a power-law of t~ 1,2 for delay times > 400 Myr. 

2 The supernova rate in the lo cal Universe as a function of Hubble 
type was recently presented in ILi et all j2010al ) who found the SN 
la rate to be constant across galaxy Hubble type, with a value of 
0.136±0.018 SNuM (1 SNuM = SNe la per century per 10 10 M ). 



Thermonuclear double- detonations as possible SNe la 3 



Ha to be detectable in the nebular spectra more frequently 
(|Leonardll2007l ; lHavden et alj|2010h . 

The DDS is an attractive model for SNe la, given the 
theoretically predicted occurrence rate as well as the fact 
that CO WD mergers are systems that are essentially devoid 
of hydrogen. The main argument against the DDS is that de- 
tailed WD merger calculations between two CO WDs with a 
total mass >Mch indicate that the merging process, while it 
can potentially lead to a th ermonuclear explo sion if the cor- 
rect conditions are satisfied (|Yoon et al.ll2007T ) , is more likely 
to result in collapse an d form a neutron s t ar; an accretion 
induc e d collapse (AIC, Mivaii et al.l Il980l ; ISaio fc Nomotol 
1 19851 : iNomoto fc Kondol 199ll . and re ferences therein) . 
Based on modern c ollapse calculations dFrver et al. 19991: 
iDessart et all 120071 : lAbdikamalov et alJl20ld ). iFrver et all 
l|2009h found that these AICs produced outbursts that were 
~l-3 magnitudes dimmer than typical SNe la, arguing that 
AICs could only explain a few abnormal la explosions. 

In a merger of two CO WDs, once the larger (less mas- 
sive) WD fills its Roche-lobe, it is likely to be disrupted 
and rapidly accreted by the companion. This process can be 
quite violent, an d might under the righ t conditions lead to a 
SN la explosion (|Piersanti et al.ll2003l ). For example, in the 
DDS case involving the merger of two WDs with a mass ra- 
tio cl ose to unity and WD masses ~0.9 Mq (|Pakmor et al.l 
I2010T ), critical conditions for the successful initiation of a 
detonation ([Seitenzahl et alj 12009] ) can be obtained. The 
iPakmor et al. I (|20ld ) study found that these DDS systems 
can both in number and in observational characteristics ac- 
count for the population of sub-luminous 1991bg-like SNe la. 
However, the 1991 bg-like systems only account for a small 
fraction of SNe la (|Li et allboiobh . 



For lower mass ratios fsee lPakmor et al]|2010l . Pakmor 
et al. 2011) it is unlikely that the merger would lead to 
a SN la as the achieved densities are not high enough to 
enable a detonation to occur. In such WD mergers, high ac- 
cretion rates onto the relatively 'cold' primary WD can lead 
to carbon burning off-centre, where the densities are too 
low, and carbon does not burn explosively. The primary CO 
WD will in turn burn carbon an d evolve into an oxygen - 
neon-magnesium (ONeMg) WD l|Nomoto fc Kondolll99ll ). 
As the ONeMg WD increases in mass, density and temper- 
ature conditions become more favourable for electron cap- 
tures, which in turn remove electron degeneracy pressure 
from the undisrupted WD. As the WD approaches Mch, the 
central densities continue to increase and the WD collapses 
to become a neutron star before a thermonuclear explosion 
can take place. 

Despite this, in population synthesis calculations it is 
typically assumed that all mergers of CO-CO WDs produce 
a SN la provided that the total mass exceeds Mch- If some 
(or many) of these CO-CO mergers result in AIC then the 
observed SN la rates cannot be fully explained by the DDS 
model. Thus, if it is true that the majority of the DDS sys- 
tems cannot produce events that look like SNe la, and there 
are not enough SDS or HeRS events, then a significant frac- 
tion of SNe la remain to be accounted for. 



2.3 The helium donor formation channels 

2.3.1 Chandrasekhar mass explosions 

We delineate between the SDS and HeRS since the latter 
can involve either one degenerate star where the donor is 
helium-burning, or two, where the donor is a helium-rich 
WD. Double WDs will have very close orbits (orbital peri- 
ods <70 minutes) as is the case for typical AM CVn bina- 
ries fsee lNelemans et alJliooH . |2010| . for a discussion on AM 
CVn stars). In Paper I the rates and delay times from three 
formation channels involving exploding WDs with masses 
^ M Ch (DDS, SDS and HeRS) were investigated. Both 
helium-donor channels were referred to as the 'AM CVn 
channel' in Paper I. In this paper we adopt the acronym 
HeRS for all SN la progenitors in which the Chandrasekhar 
mass WD explodes once it has accreted sufficient mass in 
stable RLOF from a helium-rich companion, whether the 
donor is degenerate or non-degenerate. This scenario in- 
cludes AM CVn binaries as well as WDs accreting from all 
helium-burning stars. 

SN la rates of the HeRS scenario leading to 
SNe la have been previous l y inv es tigated by differen t 
groups: ISolheim fc Yungelsonl j|2005l) ; iRuiter et alj (|2009h ; 
IWang et all (|2009al lbf): lMeng fc Yangj (|2010l ). some of whom 
considered only the helium-burning star channel. In the ma- 
jority of studies it was found that the HeRS is unable to ac- 
count for the rates of SNe Ia0 In many cases, theoretically- 
motivated studies of the HeRS channel produce SNe la with 
short delay times (< a few hundred Myr), and are not able 
to account for a large number of systems at long delay times. 



2.3.2 sub -Chandrasekhar mass explosions 

Thermonuclear explosions may occur in systems with a 
sub-Mch (probably CO) WD accr eting via stable RLOF 
from a helium-rich companion (llben fc Tutukovl Il99ll : 
iTutukov fc Yungelsonl 19961 ; lYungelson fc Livioll2000l ). It has 
been shown that at certain (low) accretion rates on to 
the WD helium flashes on the WD sur face are inhibited 
|Kawai et al.ll 19871 ; Hvanova fc Taamll2004 ). and the WD can 
steadily and efficiently build up a massive layer of helium on 
the surface. In such a massive de generate helium shell (e.g., 
~0.1 M Q of helium. iTaaml 1 1980T ). a flash may likely evolve 
as a violent detonation, which may also trigger a detonation 



complete star (e.g., 


Livnelll990l; IWooslev & Weaver! Il994l; 


iLivne & Arnettj|l995 


) 4 Tutukov & Yungelson ( 19961) found 



3 IWang et al. | ll2009ah found Galactic SN la rates which were 
higher than the other studies: ~10~ 3 yr -1 . However this rate is 
likely somewhat optimistic because they consider a rather large 
range of orbital periods at the moment of RLOF onset between 
the WD and the helium star (up to >100 days for the most mas- 
sive WD accretors). 

4 It is possible that the companion can be hydrogen-rich where 
hydrogen burns steadily on the surface of the WD, building up 
a helium-rich l ayer on top of the CO WD which can detonate 
jKenvon et alj 1 19931 : iPiersanti et alj Il999l : lYungelson k. Livid 
l2000f) Such a double-deto nation scenario for h ydrogen-rich 
donors was investigated by lYungelson et al" I Jl995h . who have 
shown that sub-Chandrasekhar mass WDs accreting hydrogen in 
symbiotic binaries may be capable of producing up to ~ i of 
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that the rate of sub-Mch SNe la from the non-degenerate 
helium star channel might be high enough to account for 
the Galactic rate if these explosions are comparable in lu- 
minositiy to normal SNe la, though these types of events 
were found to have somewhat short delay times and cannot 
account for the number of SNe la in old stellar populations. 

Sub-Mch models of SNe la are appealing for a number 
of reasons. Population synthesis calculations have already 
shown that the DTD for the double-detonation sub-Mch ex- 
plosion m odel (sometimes referred to as edge-lit detonation, 
ELD; see lYungelson fc Liviol |2000|. and references therein ) 
spans a wide range at early times. Yungelson fc Livicl (|2000l ) 
investigated the double-detonation scenario for helium star 
and hydrogen-rich donors, and found a corresponding de- 
lay time of ~ 30 Myr— 1.3 Gyr for the helium star donor 
systems. Also, population synthesis modelling has indicated 
that the number of potential progenitors for the double- 
detonation sub-Mch SNe la involving helium-rich donors 
alone may be large enough to account for the observed rate s 
(e.g jTout et alj200ll ; lRegos et alj2003l : lRuiter et alj|2009l ). 
Additionally, simplistic studies of pure detonations of sub- 
A/ch WDs (in the absence of any overlying helium shell) in- 
dicate that the synthetic light curves and spectra from sub- 
Mch explosions may be able to reproduce a surprising num - 
ber of the observed properties of SNe la (|Sim et al.ll2010l ). 
and even variations within the class that could be associated 
with diff erences in t he ma ss of the exploding WD. 

The lSim et ail l|2010l ) work, however, neglects the issue 
of how the outer helium layer will affect the observables. Sev- 
eral previous studies have calculated detailed synthetic light- 
curv es and spectra of sub-Mch double-detonation models 
(e.g. iHoeflich fc Khokhlwll 19961 ; iNugent et al.lll997l ). These 
concluded that such explosions would likely not lead to 
events with observational properties characteristic of normal 
SNe la. In general, the light curves were found to rise and fall 
too rapidly compared to 'normal' SNe la while their spectra 
were too blue to match sub-luminous SNe la and lacked suf- 
ficiently strong features of intermediate mass elements, such 
as Si and S. Importantly, most of these discrepancies with 
observation can be traced to the presence of the products of 
helium burning ( 56 Ni and other iron-group elements) in the 
outer regions of the ejecta, and those studies mainly consid- 
ered systems in which a relatively massive (~0.2 Mq) he- 
lium layer had accumulated on the W D (~ 0.6 Mq) surface. 
More recently, iBildsten et all l|2007l ) and IShen fc Bildstenl 
l|2009D have shown that conditions suitable for detonation 
in the WD might be reached for somewhat lower helium 
shell masses than considered in most previous studies: per- 
haps as low as 0.05 Mq for a CO WD (core) mass of 1.0 Mq 
(in general the more massive the CO WD, the less massive 
the accumulated helium layer needs to be for a detonation). 
iFink et all (j201fj| ) have shown that, even for such low helium 
shell masses, detonation of the helium will robustly lead to 
an explosion of the underlying WD. With a significantly 
lower He shell mass (and thus fewer iron-group elements 
in the outer ejecta), this may open the door for double- 



all SNe la, provided that accreting WDs with masses as low as 
0.6 Mq are able to successfully undergo double-detonations. How- 
ever in this work, we only consider helium-rich donors as possible 
companions for double-detonation sub-Mch SN la progenitors. 



detonation sub-Mch models whose spectra and light curves 
are in better agreem ent with obse r ved S Ne la. This has 
been investigated by iKromer et all (|2010|) w ho computed 
synthetic observables for the lFink et al. nfioiop simulations. 
They showed that even very low mass (0.05 Mq) helium 
shells affect the observable display and can lead to spec- 
troscopic signatu res that are not char acteristic of observed 
SNe la. However. I Kromer et al.1 (|2010l ) also highlighted that 
the results are highly sensitive to the details of the nucle- 
osynthesis that occur during burning of the helium shell. 
Modifications to the burning - as might be achieved by con- 
sidering a composition other than pure helium - could al- 
low the model predictions to achieve much better agreement 
with observation. 

Taken together, the body of theoretical work strongly 
suggests that the sub-Mch double-detonation scenario is 
physically realistic. Depending on the details of the accumu- 
lated helium layer and its burning products, the explosion 
may closely resemble observed "normal" SNe la or it might 
be highly spectroscopically peculiar - but regardless of this, 
it certainly can be bright enough to be readily observable - 
for a CO WD of around 1.0 Mq the luminosity produced fol- 
lowing detonation is expected t o be close to tha t of a normal 
SN la (|Shigevama et al.1 1 19921 ; ISim et al.ll2010l ). Given that 
potential progenitors are also expected to be common, we 
are therefore compelled to further investigate this progenitor 
scenario. If these explosions can produce events that resem- 
ble "normal" SNe la then it is of interest to quantify the 
fraction of observed SNe la that might be accounted for via 
this channel. Alternatively, if these explosions are realised 
in nature but are spectroscopically peculiar, it is important 
to estimate their predicted rate and consider whether the 
apparent lack of observational detections is a major concern 
for the established theory; the lack of such events may chal- 
lenge our understanding of either (or both) the explosion 
physics or the progenitor binary evolution. For this scenario 
it is of particular interest to consider the DTD predicted for 
this class of explosion and to investigate any correlations be- 
tween the properties of the exploding system (particularly 
the mass of the primary WD, which likely determines the 
brightness of the explosion) and the age of the stellar pop- 
ulation in which it resides. 



3 MODEL 

It has been shown that a WD accumulating helium-rich 
material may be capable of exploding as a SN la if the 
correct co nditions are satisfied, even if the WD is be- 
low M r. h jTaaml Il980l; llben et all Il987l; llben fc Tutukovl 



199ll: IWooslev fc Weaver! 19941; llnvne fc AmettJ Il995l : 
Ivanova fc Taaml l2004l )"~In iBelczvnski et al.1 (|2005l ). sub- 



Mch models were calculated, though DTDs for different for- 
mation channels were not discussed separately, and rates 
were not presented in that work. In Paper I, rates and delay 
times were presented only for >Mcjh WD mass models where 
as here, we additionally include sub-Mch SNe la progenitors 
in our study. 

All sub-Mrjh SN progenitors in our calculations involve 
a CO WD accreting via RLOF from a helium-rich compan- 
ion. As in Paper I, if the donor is a WD then it can be either 
a helium WD or a hybrid WD; a WD with a CO core and 
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a helium-rich mantle (e.g.. iTutukov fc Yungelsonlll996l . and 
references therein). Hybrid WDs are formed through binary 
evolution when a red giant is stripped of its envelope through 
binary interactions. In cases where the stripped helium core 
does not reach the helium-burning phase, a helium WD is 
formed. 

In the following sections we compute and discuss rates 
and delay times for the aforementioned SN la evolutionary 
models that have been proposed as the most promising for- 
mation channels for SNc la: 

• the DDS 

• the SDS 

• the HeRS 

• the double-detonation sub-Mch scenario involving 
helium-rich donors 

It is important to keep in mind that, while evolution 
of close binaries remains an active field of research and dis- 
covery, no concrete constraints currently exist for the evolu- 
tion of mass-transferring binaries, nor for the common en- 
velope (CE) phase. The CE phase is certainly one of the 
most poorly understood phenomena in close binary evolu- 
tion, and a theoretical picture of CE evolution is not yet 
available. Since there are a limited (though growing) number 
of observations available to guide our choice of parameters, 
we present results for three different common envelope re- 
alizations which most effectively bracket the uncertainties. 
For the growth of CO WD s during stable RLOF, as was 
done in lRuiter et alJ (j2009h we present the results from our 
population synthesis model using a detailed WD accretion 
scheme, which was constructed by adopting various input 
physics from the literature. 

We use t he StarTrack populatio n synthesis binary evo- 
lution code (|Belczvnski et af] l200ct ) to evolve our stellar 
populations employing Monte Carlo methods. The code 
has undergone m any revisions since the first code descrip- 
tion publication l|Belczvnski et alj|200fl ). Many of the up- 
dates concerning accretion on to WDs can be found in 
iBelczvnski et alj (|2005l .[ 20081 ) , though since then we have in- 
corporated an updated prescription for accr etion of hydro- 
gen o n WDs by incl uding calculations f r om Nomoto et all 
(20071 ) in addition to lPrialnik fc Kovetj (|l995h . The initial 
distributions for binary orbital parameters (orbital periods, 
mass ratios, etc.) are the same as described in Paper I, sec- 
tion 2. 

In Paper I, it was assumed that the ejection of the enve- 
lope of the mass-losing star during a CE phase came at the 
expense of removing the orbital energy of the binary, as dic- 
tated by t he well-known ' energy-balance' (or 'a-formalism') 
equation (|Webbink|[T98l ). with qce representing the effi- 
ciency with which the binary orbital energy can unbind the 
CE, and A is a para metrization of the structure of the donor 
star (|de Koollll990h ; both qce and A are fairly uncertain. 

For Models 1 and 2 from Paper I, «ce x A values of 
1 and 0.5 were adopted, respectively. The major difference 
was that Model 1 (more efficient removal of the CE) resulted 
in an overall higher number of SNe. In the current paper, we 
keep all model parameters the same as in Model 1 of Paper I 
for one model; we refer to this model as Model Al (standard 
model). However, in order to explore the sensitivity of the 
physical mechanism of CE ejection, which is still not un- 
derstood, we have run two additional sets of models. There 



has been some recent observational |Zorotovic et alj|20ich 
as well as theoretical |Passev et al. I l20ld . in prep.) evidence 
that the value for qce lies between 0.2 and 0.3. Addition- 
ally, for low-mass stars a v alue of 0.5 is often adopted for A 
(van der Sluvs et"afll2006h . Thus to best bracket our uncer- 
tainties for the energy balance prescription of CE evolution, 
in a second model we employ a very low CE ejection effi- 
ciency: qce = 0.25 and A = 0.5 yielding qce x A = 0.125. 
We refer to this model as Model A. 125. In a third model, we 
assume a different parametrization for the treatment of the 
CE phase. We employ the '7' prescription for CE evolution 
|Nelemans et alj|200d ) every time a CE event is encountered 
in the code. Although the combination '07' prescription is 
considered to be the preferred prescription by some groups, 
the two formalisms stress different physics (one focusing on 
energy conservation, one on angular momentum balance), 
and we do not mix the formalisms in the current paper. 
Henceforth we refer to the 7 model as Model G1.5. For 
Model G1.5, all physical parameters are identical to Mod- 
els Al and A. 125, except that when unstable mass transfer 
is encountered and a CE ensues, the orbital separation of 
the binary changes not as a consequence of removing grav- 
itational binding energy from the orbit, but linearly as a 
function of mass loss (and hence angular momentum loss), 
parametrized by the factor 7: 



Of 



l- 7 



Mo 
Mto< 



2 . \ 2 

Mtot.f / Mdon.i Mcom 



Mu 



Mdon.f Mo, 



(1) 



where Mdon.i is the initial mass of the (giant) donor star 
just prior to the CE, M e j is the ejected mass (assumed to 
be the mass of the giant's envelope), M CO m is the mass of 
the companion (assumed to be unchanged during the CE), 
Mdon.f is the final mass of the donor once the envelope has 
been ejected, <Zi is the initial orbital separation, af is the 
final orbital separation, and M to t,i and M to t,f represent the 
total mas s of the binary before and after CE, respectively. 
Following iNelemans et all (|200Ch . we have chosen 7 = 1.50 



3.1 Sub-Chandrasekhar mass model: Assumptions 

We adopt the prescription of llvanova fc Taaml (|2004 ). ap- 
plied to accretion from helium-rich companions only, to de- 
ter mine when a particular binary undergoes a sub-Mch SN 
la l|Belczvnski et all 120081 . see section 5.7.2 for equations). 
In short, we consider three different accretion rate regimes 
for accumulation of helium-rich material on all CO WDs, 
a dopting the input physics for helium accretion on to WDs 
of lKato fc Hachisul (|1999l . l2004h . Clearly, our results are sen- 
sitive to the adopted helium accumulation efficiency. It may 
also be possible that for the less-common hybrid donors, 
the star's helium-rich envelope becomes stripped in RLOF 
before a helium shell detonation on the CO WD occurs. 
This would have an effect on the purity of the helium-rich 
shell, though we do not investigate the consequences here. 
We assume that mass is exchanged conservatively in all cases 
where the transferred material would be CO-rich. 



5 We note here that the 7 CE equation in[Bclczyiiski et al.l ||2008| . 
equation 55) is missing an exponent, though the CE evolution 
is properly carried out in the StarTrack code with the correct 
equation. 
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At high accretion rates (~ 10~ 6 M@yr , and for all 
cases where the WD accretor mass is < 0.7 M©), helium 
burning is stable and thus mass accumulation on the WD is 
fully efficient (r/ acu = f). At somewhat lower accretion rates 
helium burning is unstable and the binary enters a helium- 
flash cycle, thus accumulation is possible but is not fully 
efficient (0 < ?7acu < 1). In both of these aforementioned 
accretion regimes, the CO WD is allowed to accrete (and 
burn) helium, and its total mass may reach Mch and explode 
as a SN la through the HeRS channel. However, for low 
accretion rates (~ 10 -8 Mgyr -1 ), compressional heating 
at the base of the accreted helium layer plays no significant 
role, and a layer of unburned helium can be a c cumu lated 
on the WD surface. Following llvanova fc Taaml l |2004 ). we 
assume that if such a CO WD accumulating helium enters 
this 'low' accretion rate regime and accumulates 0.1 Mq of 
helium on its surface, a detonation is initiated at the base 
of the helium shell layer. Consequently, a detonation in the 
core of the CO WD is presumed to follow, and we assume 
that a sub-Mch SN la takes place. Only accreting WDs with 
a total mass ^0.9 Mq are considered to lead to potential 
sub-Mch SNe la in this work, since lower mass cores may 
not detonate, and if they do they are unlikely to produce 
enough radio active nickel and hence will not be visible as 
SNe la (e.g., ISim et al.l [20101 . table 1). Thus in all future 
discussions we refer to sub-Mch systems whose total WD 
mass (CO core + helium shell) is at least 0.9 Mq at the time 
of SN la unless otherwise noted; for our population synthesis 
model, this intrinsically implies that all exploding sub-Mch 
SNe la have CO WD 'core' masses >0.8 Mq. Helium-rich 
WDs are simply not massive enough, and we assume that 
ONeMg WDs do not make SNe la. 



4 RESULTS 

Here we present the DTD and rates for all of our SN la mod- 
els, as well as CO WD core masses for our sub-Mch models. 
Our results are discussed in the following subsections, but 
here we give a brief outline of our findings: 

We have investigated the DDS, SDS, HeRS, and the 
sub-Met scenario for three different CE realizations. Within 
the framework of our adopted models, we find that only two 
SN la formation scenarios are capable of matching the ob- 
served SNe la rates: the DDS and the sub-Mch channels. 
The most favourable model in terms of matching observa- 
tional rates is model Al (a X A = 1). For models Al and 
G1.5 (7 = 1.5), the adopted sub-Mch scenario is dominant 
at nearly all epochs <5 Gyr, however the sub-Mch chan- 
nel rate is too low for our low-efficiency CE model (A. 125). 
For Model A. 125, no single progenitor, nor an admixture of 
all of the progenitors combined are able to account for the 
observed rates of SNe la. 



4.1 Delay times 

In Figure Q] we show the delay time distribution of the 
four aforementioned progenitor channels for Models Al, 
A. 125 and Model G1.5. We note that the bumpiness in the 
smoothed plot is due to Monte Carlo noise. For our DTD 
normalisation of all models, we have assumed a binary frac- 
tion across the entire initial stellar mass function of 50% (| 



of stars are in binaries), and we show the DTD normalised 
to stellar mass (SNuM and SNe yr~ 1 Mg 1 ). The mass rep- 
resents the mass in formed stars, which includes mass which 
has potentially been expelled from stars in SNe or thermal 
pulses for example. In section 4.2, we give the delay times 
in tabular form. 

Along with our theoretical DT Ds, we show the obs erved 
(cosmic) DTD from the literature (|Maoz et alj|2010a[ ). We 
wish only to compare the relative DTD shapes and not the 
absolute rates, since the normalisation of our StarTrack 
DTD differs substantially from the normalisation techniques 
used in recovering the various observational DTDs. The dif- 
ference between the (higher) rates of observed SNe la and 
the rates from population synthesis is visible, and the ap- 
parent discrepancy is not yet fully resolved. It has been 
suggested that binary population synthesis codes tend to 
under-predict the SN la rates compared to the rates infered 
from recent observations, though one must keep in mind 
that many uncertainties are associated with the DTD recov- 
ery methods, i.e., extinction, star formation history, and the 
use of spectral populatio n synthesis codes which neg l ect th e 
exis tence of binaries (see |Pe Ponder fc Vanbeverenl l|2004 ). 
also lEldridge fc Stanwavl (|2009h have found that inclusion 



of massive binaries in spectral synthesis codes plays an im- 
portant role in recovering accurate host galaxy properties). 

Model Al. As was found in Paper I Model 1, the DDS 
distribution for Model Al (top panel of Figure [T]) follows a 
power- law distribution with (see also Figure [2]), while 

the SDS distribution is somewhat flat with no events with 
delay times less than 460 Myr. The reason why the SDS 
does not harbour very prompt events is directly linked to 
the donor star's initial ZAMS mass. When the secondary 
ZAMS mass is >2.8 Mq, the binary will enter a CE phase 
when the secondary fills its Roche- Lobe, rather than a stable 
RLOF phase. In such a case, the binary will not become an 
SDS SN la, though may under the right circumstances evolve 
to SN la from the HeRS channel. The SDS events at long 
delay times originate from progenitors with very low-mass 
MS donors, which take many Gyr to evolve to contact under 
the influence of magnetic braking. The HeRS DTD consists 
mostly of systems with relatively short (~100 Myr-2 Gyr) 
delay times, with very few events at longer delay times. We 
refer the reader to Paper I for a description of these DTDs. 

The sub-Mch systems can easily by eye be grouped into 
two classes: those prompt SNe which occur with delay times 
<500 Myr, and those with delay times above ~800 Myr, with 
very little overlap. Not surprisingly, these two classes of SNe 
la stem from two very different formation scenarios. Those 
with short delay times consist of progenitors which involve a 
helium-burning star donor, whereas the rest mainly consists 
of helium WD donors (systems with hybrid WD donors span 
~0.3-3 Gyr delay times). We find that progenitors with he- 
lium star, helium WD and hybrid WD donors comprise 13%, 
78% and 9% of SNe la, respectively. We note that ~35% of 
sub-Mch SNe la explode within 1 Gyr of star formation. 

The prompt component accounts for 13% of all sub- 
Mch SNe la that explode within 13 Gyr of star formation. 
Nearly all of these systems (96%) have helium star donors, 
with the rest having hybrid WD donors. The delay time is 
governed by the MS lifetime of the donor star. The com- 
panions with ZAMS masses >3 Mq evolve off of the MS 
within <400 Myr. After the first CE, which leaves behind a 
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Figure 1. Lines represent the DTD for SNe la. Top panel: Model 
Al. Middle panel: Model A. 125. Bottom panel: Model G1.5. The 
number of SNe la per year per unit stellar mass born in stars (at 
starburst t = 0, 50% binarity) is shown for the DDS (blue), SDS 
(red), HeRS (green), and sub-Mch (magenta) channels. The sub- 
Mch SN la DTD clearly shows two distinct populations for Mod- 
els Al and G1.5: the helium star channel (spike at delay times 
<500 Myr) and the WD channel (from ~800 Myr to a Hubble 
time). The helium star channel however is absent in Model A. 125. 
In the top pane l we ad ditionally show the DTD [SNuM] compiled 
by iMaoz et al. | i2010al , table 1), which is fit relatively well by a 



power-law ~ t 



The data points showing the observed DTD 



are computed using a different normalisation technique (see text), 
and thus we show the points for comparison of the DTD shapes 
and not the absolute numbers. We note that assuming a different 
binary fraction or IMF would change the level of our normalisa- 
tion. 



CO primary WD and a MS secondary star, the secondary 
(e.g., on the Hertzsprung gap) will fill its Roche-Lobe and 
mass transfer is once again unstable leading to a second CE 
phase. The CE leaves the CO WD and newly-formed naked 
helium star on a close orbit (^35-40 min). Within a few 
Myr, the orbit decreases to ~25 min, and the helium star fills 
its Roche-Lobe. However, initial mass transfer rates for the 
helium star channel are low enough to enable accumulation 
of the helium shell to commence immediately: typically such 
systems have initial mass transfer rates ~2 x 10 -8 Mq yr _1 
(for a discussion on the evolution of low mass helium stars 
in accreting binaries see lYungelsonlkoOSl ). 

The delayed component (delay times >500 Myr) com- 
prise the other ~87% of the sub-Mch progenitors. Binaries 
with helium WD donors make up 90% of the delayed com- 
ponent, while 10% have hybrid donors. These binaries also 




log(delay time/Myr) 



Figure 2. Delay time distribution for the DDS (blue) and sub- 
Mq, (magenta) channels for model Al (standard). We show two 
power-laws alongside the DTDs: the DDS is relatively well-fit by 
a power-law t , where as the sub-Mch model closely follows a 
power-law with t~ 2 beyond 1 Gyr, where all progenitors have 
helium-rich WD donors. 



evolve through two CE phases, as is expected for the evolu- 
tion of AM CVn binaries. Similar to the DDS, the time- 
scale governing the DTD for the helium WD channel is 
largely set by the gravitatio nal radiation time-scale (see also 
iTutukov fc Yungelsonll 19961 ). However unlike the DDS, these 
WDs do not merge upon contact, but enter a stable phase 
of RLOF. Like the DDS DTD, the sub-M Ch DTD follows a 
power-law above 1 Gyr, however with a steeper functional 
form of t~ 2 (Figure [21 see also section 5.2). 

Model A. 125. In the middle panel of Figure [T] we 
show the DTD for Model A. 125. Contrary to the standard 
model, the DDS DTD is lacking progenitors at longer de- 
lay times since on average the time a progenitor spends as 
a detached double WD is decreased in this model (smaller 
orbital separation following the CE phase). For this model 
the SDS progenitors can have shorter delay times compared 
to the standard model due to the fact that the post-CE 
separations are overall smaller. Thus, a low CE efficiency 
model is more favourable for the production of SDS SNe la. 
The HeRS channel has some very prompt events (helium 
star channel), although the low CE efficiency serves to re- 
sult in a merger during CE more frequently than in Model 
Al. The events at delay times ~1~2 Gyr belong to progen- 
itors with helium WDs, while events at long delay times 
(> a few Gyr) also involve helium WDs but belong to the 
evolved low-mass MS donor channel discussed previously. 
The DTD of the sub-Mch progenitors looks drastically dif- 
ferent from that of the standard model, and lacks a prompt 
component. This model is the only of the three which does 
not display a prominent division of the sub-Mch progenitor 
channels; in fact there are no sub-Mch SNe la originating 
from the helium star channel, since those progenitors will 
encounter unstable RLOF too early in their evolution. With 
the adopted CE prescription, it is very difficult (or impossi- 
ble) to produce helium star donor channel sub-Mch SNe la 
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within our model framework, and thus there are no prompt 
events (the first SN la from the sub-Mch channel occurs at 
-1.7 Gyr). 

Model G1.5. In Figure [1] bottom panel, we show the 
DTD for Model G1.5. Gravitational radiation plays a less 
significant role for the DDS since following the CE phase 
the binary orbit is still rather wide. Similar to the other 
two models, the DDS contributes the majority of its events 
at very early times followed by a decline. The SDS channel 
displays no prompt events, because the first CE event does 
not lead to a dramatic decrease in orbital separation. In 
general, the SNe la with short SDS delay times from Model 
Al will evolve into detached double CO WDs in Model G1.5, 
since the binary orbit is not small enough for mass transfer 
to begin once the secondary evolves off of the MS and fills 
its Roche-Lobe. The HeRS channel leads to SNe la with 
very short delay times, though there are events at long delay 
times but their frequency is for the most part too low to be 
seen on the figure. The sub-Mch DTD has the same general 
shape as Model Al: the prompt and the delayed components. 
We note however that SNe la with delay times less than 
1 Gyr follow a different evolutionary sequence compared to 
the corresponding events of Model Al. In Model Al, the first 
mass exchange interaction occurs when the primary is an 
AGB star, where as for Model G1.5 the first mass exchange 
event (CE or stable RLOF) occurs when the primary is less 
evolved; a sub-giant or giant. This occurs since the semi- 
latera recta (and thus in general, the separations) of the 
G1.5 sub-Mch progenitors which explode in our models are 
smaller when mass transfer begins, as well as the fact that 
the primaries for this model are somewhat more massive 
than compared to the standard model and thus they evolve 
more quickly 

4.2 Rates 

In Table 1, we show the DTDs in tabular form (rates as a 
function of epoch) for our models. We estimate the Galactic 
SN la rate by convolving the DTD (in units of SNe/time) 
with a constant star formation history from — 10 Gyr with 
a total mass born in stars of 6 x 1O 1O M0; see section 4 of 
Paper I. In this paper, we do not quantify the Galactic rate 
estimates explicitly since imposing a particular star forma- 
tion history serves to add sources of uncertainty to our DTD 
calculation, however we give some numbers as a guide. In 
principle, all of the important information is already pre- 
sented in the DTD plots and Table 1: it is possible to con- 
volve the specific DTD with any star formation history of 
choice in order to achieve a particular SN la rate for a given 
stellar population. 

Model Al: This model produces the highest number of 
SNe la out of our 3 models. Table 1 (left) is very similar to 
table 1 (Elliptical column) in Paper I, though here for all 
tables we additionally include the rates of sub-Mch SNe la, 
as well as two additional epochs: 0.1 and 1 Gyr after star for- 
mation. Slight variations between the numbers in this study 
and table 1 of Paper I are due to a slight increase in vol- 
ume of data, and thus a reduction in noise from low-number 
statistics. We find that the rate of our adopted sub-Mch 
SN la model exceeds all other progenitor channels between 
~0.7 and 5 Gyr, and these systems are enough to account 
for the observed SN la rate, with a calculated Galactic rate 



Table 1. Rates of SNe la (SNuM, 50% binarity) for the four 
progenitor formation scenarios considered in this work, following 
a starburst at t = 0. Models Al (left), A. 125 (middle) and G1.5 
(right). 
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of ~2.6 x 10~ 3 SN la yr _1 (including all systems with a to- 
tal WD mass >0.9 Mq). For comparison, the DDS rate is 
—2 x 10 -3 SN la yr -1 . Both of these values are within the 
estimate from ICappellaro et al] (jl999h of 4 ± 2 x 10~ 3 SN 
la yr _1 . As was determined in Paper I, the Model Al DDS 
rates are able to (just) account for the observed Galactic 
rate of SNe la, whereas both the SDS and HeRS channels 
fall short by over an order of magnitude. 

Model A.125: This model produces the least SNe la pro- 
genitors out of our three models. The DDS is significantly 
decreased in number (Table 1, middle), but is still the dom- 
inant channel at most times under a few Gyr. The SDS rate 
exceeds the DDS rate above 3 Gyr, though the overall rates 
are still too low for any progenitor in this model to account 
the observed SN la rates. The rates of the HeRS SNe la are 
too low; many binaries do not survive both CE events to 
become progenitors. Similarly, the sub-Mch progenitors are 
not easily formed in this Model. 

Model G1.5: The overall rates for this model (Table 1, 
right) are lower than found in the standard model, though 
not as low as found for Model A.125. In the DDS, since 
the binaries take a longer time to reach contact (e.g., it can 
easily be more than a Hubble time), the overall SN la rates 
are rather low compared to Model Al with an estimated 
Galactic rate of —2 x 10~ 4 yr _1 , which is about a factor 
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of 10 too low. The SDS channel produces very few events 
before 2 Gyr, and matches those of the DDS at ~3 Gyr, 
while the HeRS channel produces events with delay times 
<1 Gyr and few events at later times. Even though the sub- 
Mch DTD exhibits the same general shape as found in the 
standard model, the rates are overall too low being roughly 
comparable to those of the DDS of this model (Galactic rate 
estimate ~3 x 10~ 4 yr _1 ). 

4.3 CO core masses 

In the sub-Mch scenario, the brightness is expected to be 
largely determined by the mass of the underlying CO WD. 
In Figure [3] we show the mass of the CO WD 'core' (to- 
tal WD mass minus the helium shell mass) at time of SN 
la. As mentioned previously, a detonation of a ~0.7 Mq 
core WD would likely not look like a normal SN la. Since 
we currently lack a theoretical lower mass limit for which 
exploding CO core masses could potentially exhibit features 
which are characteristic of SNe la, for completeness we show 
the CO core mass at explosion for the entire mass spec- 
trum for exploding sub-Mch cores. We draw a vertical line 
at M core = 0.8 Mq, above which the systems are considered 
to be sub-Mch SNe la in our models. 

The core mass distributions look very different for all 
three models. In the top panel of Figure O we show the core 
mass distribution for Model Al. The progenitors of bina- 
ries with low core masses (< 0.7 Mo) go through a differ- 
ent evolutionary channel than those with higher core masses 
since they start out with smaller semi-latera recta and only 
evolve through one CE event. The cores associated with the 
helium star channel span both low and high masses, though 
for our adopted sub-Mch scenario they have slightly higher 
core masses on average compared to the WD channels. The 
hybrid WD channel shows a similarly flat distribution, which 
is not unexpected since many of these systems undergo an 
evolutionary sequence which is like that of a typical pro- 
genitor from the helium star channel. For the helium WD 
channel which comprises the majority, the masses decrease 
fairly steadily in number with increasing mass, since there 
are simply a larger number of less-massive CO WD cores to 
start with. There is a clear lack of CO core masses below 
~ 0.7Mq. Typically these CO core progenitors will accrete 
(and burn) at least 0.1 Mq (often ~ 0.2 Mq) of helium at a 
high accretion rate before the phase of helium accumulation 
begins for the 0.1 Mq shell, and thus we find no CO cores 
from this channel with very low masses. However, there are 
a number of exploding cores with masses ~ 0.7 — 0.8 Mq. 
One has to also consider the possibility that a low-mass 
(< 0.8 Mq) CO core + helium shell may not reach suffi- 
cient conditions for a detonation to take place, which might 
explain why we would not observe a large number of these 
events. 

For Model A. 125 (Figure [3j middle panel), the separa- 
tion between the helium star and double WD channels is 
quite distinct. The lowest mass CO cores belong to progen- 
itors with helium star donors, and in our adopted sub-Mch 
model all of these binaries have CO core masses which are 
too low to qualify as SNe la. The ZAMS masses of these CO 
cores are small, ~1. 8-2.1 Mq, and these stars are unable to 
build a massive CO core before the first CE is encountered. 
Additionally, binaries which start their final RLOF phase 
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Figure 3. Distribution of Star Track CO WD core masses which 
managed to accumulate a 0.1 Mq shell of helium. A double- 
detonation was assumed to follow in all cases. A vertical line 
is drawn at MfjOcore = 0.8 Mq, above which all systems 
are assumed to lead to sub-Mch SN la in our three models 
(Afcocoro + Aijicshcll = total WD mass). Top panel: Model Al. 
Middle panel: Model A. 125. Bottom panel: Model G1.5. The he- 
lium star channel is outlined in blue, the He- WD channel is out- 
lined in red and the hybrid WD channel is outlined in black. Note 
different scales on the j/-axes. 
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when the secondary is a helium star have lower initial ac- 
cretion rates (~10 -8 Mq yr" 1 ), which allows the CO core 
to immediately accumulate (not burn) a shell of helium and 
produce a SN la without the CO WD having to grow in 
mass by an extra ~0.1-0.2 Mq. 

The distribution of CO core mass for the helium WD 
channel of Model G1.5 (Figure[3] bottom panel) is very simi- 
lar to that of the standard model. However the different evo- 
lutionary sequences allowed in this model enable the forma- 
tion of more progenitors involving hybrid WD donors. The 
mass distributions for the helium star and hybrid WD chan- 
nels peak between 0.85 and 0.9 Mq (total WD mass 0.95- 
1 Mq), which is a noteworthy feature, especially if these 
systems are shown to co ntribute to the p opulation of SNe 
la of 'normal' brightness |Sim et al.ll2010T ). 



5 DISCUSSION 

Recen t hydrodynamic explosion simulations of sub-Mch CO 
WDs (|Fink et al.ll2010T ) coupled with detailed nucleosynthe - 
sis and radiative transfer modelling jKromer fc Simll2009T ) 
have revealed that sub-Mch mass SN la models exhibit fea- 
tures whic h are characteristically similar to tho se observed 
in SNe la (|Sim et alj|20ld : iKromer et alj|20ld ). Motivated 
by these new findings, as well as population synthesis rate 
estimates, we have investigated sub-Mch SN la formation 
channels and have calculated and presented the delay time 
distribution and rates of their progenitors for three different 
parametrizations of the common envelope phase. 

We find that only the sub-Mch progenitor channel is 
able to simultaneously 

• reproduce the observed rates for our standard model 

• provide an elegant explanation for the variety among 
SN la light-curves (mass of exploding WD) 

• naturally provide a system which is devoid of hydrogen 

• produce a DTD with distinct prompt (<500 Myr) 
and delayed (>500 Myr) components, originating from 

two channels with very different evolutionary time-scales 

We think that this last point is one of the most inter- 
esting, considering the recent observational studies by dif- 
ferent groups who have found evidence for s uch a DTD 
jBrandt et all l20ld ; iMaoz fc Badenes! l20ld ; iMaoz et all 
l2010bh . 



5.1 Double white dwarf mergers: implications 

We note both the works of iGuillochon et all (|2010h , who 
investigated detonations in sub-Mch CO WDs undergoing 
rapid accretion during dynamically un stable mass transfer 
from a helium-rich WD companion, and lvan Kerkwiik et al.l 
(|20ld ). who also considered mergers of WDs with a to- 
tal mass below Mch as possible progenitors of SNe la. In 
this study we do not investigate sub-Mch mergers in de- 
tail though we briefly comment on them here. We find that 
the number of sub-Mch WD mergers in our standard model 
(considering all mergers where at least one WD is CO-rich, 
the other being CO and/or helium-rich) is nearly twice that 
of DDS mergers. While it is generally believed that a WD 
merger with a total mass below the Chandrasekhar mass 
limit would not lead to a SN la explosion, these mergers 



should produce other inter esting objects; R Coronae Bore- 
alis stars are one ex ample dWebbin k 1984; I ben et al] 199d ; 
IClavton et al. I l2007l . see also IBogomazov fc Tutukovl {{200$ )) 
and these types of merger events may be visible in upcoming 
tr ansient surveys. If we mak e a constraint similar to that 
of Ivan Kerkwiik et al.l |20ld ) counting both sub-Afch and 
super- Mch WD mergers between CO WDs with near-equal 
masses, we find that the number of mergers drops to ~42% 
of our standard model DDS rate, which is slightly too low 
to explain all SNe la. 

In our models we have assumed the commonly-adopted 
initial binary orbital configurations for population synthesis 
studies: i.e., initial separation flat in the logarithm (more 
binaries born on closer orbits relative to large orbits), and 
thus the ZAMS distribution of all semi-latera recta are the 
same for all three CE models. However, we find that for the 
low-CE efficiency case (Model A. 125), DDS SNe la progen- 
itors are only formed from systems with initial (ZAMS) or- 
bital configurations which have rather large semi-latera recta 
compared to those for our standard model. In Model A. 125, 
systems which would have made DDS SNe la in Model Al 
merge too ear ly, and never make d ouble WDs. It was already 
mentioned in iHurlev et ail |2002l ) that the initial distribu- 
tion of orbital separations in population synthesis studies 
should be distributed according to the (observed) distribu- 
tion of semi-latera recta rather than semi-major axes or or- 
bital periods alone. We note here that an initial distribution 
geared toward higher semi-latera recta than is canonically 
assumed would serve to augment the number of progenitors 
in models with low CE efficiency, making those DDS rates 
closer to those of observations. 

While the predicted rates of the DDS for our models 
do not conflict with observations, these systems are theo- 
retically expected to produce neutron stars via AIC. If this 
were the case, the AIC rate from the AlC-merger channel 
alone would be ~10 -3 per year for the Galaxy. We find 
the StarTrack AIC rate from the 'RLOF-AIC channel is 
a factor of 10 to 100 less: no more than 10 -4 per year for 
our standard model. This rate is in agreement with the up- 
per limit estimate derived from solar system abundances of 
neutron- rich isotopes, wh i ch are expected to be produced 
in AICs |Frver et al.lll999l : lMetzger et aLlbOOd ). However, if 

i) population synthesis estimates for the number of merg- 
ing CO+CO WDs with a mass above Mch are correct and 

ii) in most environments these mergers preferentially pro- 
duce AICs and not SNe la, then this could potentially be in 
conflict with the predicted abundance of neutron-rich iso- 
topes in the solar neighbourhood. However, modelling of 
AIC events, be it the 'RLOF' or 'merger' c ase, is still in 
its infancy , and many uncerta inties re main jDessart et al.l 
l200d . l2007l : lMetzger et al.ll2009l . see also lDarbha et al.ll20ld ) 
If one can say for certain that AIC events formed from the 
merger of CO WDs produce very neutron-rich ejecta, then 
this provides a potentially strong constraint on the outcome 
of these mergers; namely that a non-negligible fraction of 
SNe la must be formed through the DDS channel. On the 
other hand, it is possible that population synthesis calcu- 
lations over-predict the number of merging CO+CO WDs, 
which would also present an interesting problem for the bi- 
nary evolution community, and may challenge the idea that 
the observed ~t _1 power-law DTD of SNe la originates from 
double WD mergers alone. 
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5.2 Further remarks on delay times 

The t~ power- law shape found in the delay time study of 
iTotani et all (|200ct ) implies that the majority of progenitors 
in elliptical-like galaxies originate from binaries for which 
the DTD is most strongly governed by the time-scale as- 
sociated with gravitational wave radiation, thus these pro- 
genitors are likely to be DDS mergers (see section 3 of 
iRuiter et al.|[2003 ). However, we point out that a change 
in the CE removal efficiency a will have an effect on the 
amount of time between the last CE phase and final con- 
tact, thus affecting the shape (and perhaps to some degree 
the relevance) of the delayed DTD component. Similar to 
the DDS, our study has shown that the sub-Mch model 
DTD (e.g., Model Al) also exhibits a power-law for delay 
times >1 Gyr. This is not surprising, since the helium WD 
sub-Afch progenitors also spend an appreciable time as de- 
tached compact stars evolving to contact solely under the 
influence of gravitational radiation. However, the DTD of 
the sub-Mch ch annel falls of f more steeply than the t _1 
power-law fit of ITotani et al] (l2008h and the t -1 - 1 - i -1,3 

matching quite well 



power-law fits of iMaoz et al 



to t~ (Figure [2}. Thus, when comparing our results to 
observationally-derived DTDs, the DDS channel matches 
more closely than the sub-Afch channel. However, a very re- 
cent study of Subaru/XMM-Newton Deep Survey (SXDS) 
SNe la indicates that the DTD may be well-fit by a power- 
law of t~ 1 ' 5 (J. Okumura, private communication 20f0). It is 
of course possible that both DDS and sub-A/ch progenitors 
contribute substantially to the SN la population, potentially 
yielding a DTD of functional form somewhere in between £ _1 
and t~ 2 above 1 Gyr, which would still be in agreement with 
the majority of recent observations. 



5.3 Sub-Chandrasekhar SNe la connection to AM 
CVn stars and .la events 

Based on the observed local spac e density estimate o f 
AM CVn binaries! performed by IRoelofs et al] \200H) . 
iBildsten et al.l (|200?f ) have calculated the occurrence rate 
of the final (explosive) helium flash from ".la" systems 
in a typical E/SO galaxy with a mass of I0 11 A/ Q to be 
(7 - 20) x ICT 5 yr _1 ; i.e., 2-6% of the SN la rate in E/SO 
galaxies. .Ia events are expected to be about one tenth as 
bright as normal SNe Ia. Since our sub-Mch progenitors 
could also potentially lead to .la-like (and not SN Ia) explo- 
sions, we think it is useful to independently estimate the oc- 
currence rate for such explosions in our standard (Af ) model 
for similar (E/SO galaxy) conditions. It was already found in 
section 4.2 that the sub-A/ch rate assuming a burst of star 
formation at t = is ~4 x I0" 3 at 10 Gyr. Thus we find 
that among old stellar populations our double-detonation 
thermonucler explosions will be roughly 30 times more fre- 



6 We would like to make the reader aware of the fact that 
population synthesis studies over-predict the number of AM 
C Vn binaries in gener al co mpared to the observational results 
of IRoelofs et al.l j2007h (e.g., iNelemans et alj|200ll ; IRuiter et al.l 
|201CJ) . There may be several factors which conspire to cause the 
apparent difference, though the various possibilities are not ex- 
plored in this study. 



quent than the estimated .Ia explosion rate of lBildsten et al.l 
J2003). 

We also note that in the study of lBildsten et all (120071 ) 
it was found that the ignition mass of the helium shell in 
.las varies as a function of the underlying CO core mass 
and the rate of accretion. However for our first investigation 
of double-detonation sub-A/ch SNe we have used a more 
simplified model in which the ignition mass is always the 
same (0.1 A/©). The consequences of this on the resulting 
SN rate are not expected to be too drastic, as the time-scale 
for the helium accretion is relatively short compared to the 
evolutionary lifetime of the progenitors. This is particularly 
true for the heli um WD donor case, w hich is the scenario 
mos t relevant for Bildsten et al.l (|2007l ). 

lYoon fc Langerl ") 20041 ) found that rotation may pose a 
problem for the initiation of a detonation in accreted he- 
lium shells. They found that the spin-up of the WD due 
to the accretion and resulting dissipation due to differential 
rotation might cause helium flashes to occur for lower shell 
masses, which may lead to inhibition of a d etonation in turn 
result ing in fewer sub-Mch SNe. However, iPiro &c Bildstenl 
(|2004l ) found that the accreted material will be brought into 
co-rotation with the WD alread y at low depths with in the 
helium shell, and so as noted in feildsten et al] (|2007f ) rota- 
tion should not play a significant role in the heating of the 
helium shell and subsequent helium-ignition. 

5.4 The link between sub-Chandrasekhar SNe Ia 
and their progenitors 

While it is useful to understand how the host galaxy en- 
vironment influences the SN ejecta/observables, it is also 
fundamentally important to find a direct physical connec- 
tion between the progenitor population and the observa- 
tional characteristics of SNe Ia. For some time it has been 
known that brighter SNe Ia occur more freq uently among 
young stellar populations (|Hamuv et al.|[l995h . Could it be 
possible that sub-Afch SNe Ia arising from the (prompt) he- 
lium star channel are brighter than those from the double 
WD channel? This may be the case particularly considering 
Model Al, where the core mass of the exploding star for the 
helium star channel is on average slightly larger than for the 
double- WD channel (see Figure |3J), and thus is likely to pro- 
duce more 56 Ni. We also note that for both Model Al and 
Model G1.5, ~70% of progenitors with delay times <1 Gyr 
have CO WD masses >1.0 M Q (CO core masses >0.9 M ), 
while this fraction is only ~45-50% for progenitors with 
delay times >3 Gyr. However there is no strong trend in 
our models which indicates that more massive WDs explode 
among younger populations. The majority of the sub-A/ch 
binaries are double WDs, and the MS lifetime (ZAMS mass) 
of the primary star does not play a dominant role in setting 
the delay time. 

Another point worth considering is that the helium star 
channel progenitors undergo two CE events on a relatively 
short time-scale compared to the time the stars spend as a 
post-MS detached binary. Thus these binary systems should 
be hotter and may be more readily detectable than their 
(colder, longer lived) double- WD counterparts. Since these 
helium star channel SNe in our models are expected to oc- 
cur a few Myr after the last CE phase, the detection of such 
an explosion will probably not be inhibited by circumstel- 
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lar matter from the companion. However, since these ex- 
plosions involving helium stars are expected to be found 
among young stellar populations, they are likely to occur in 
regions of active star formation where their detection may be 
thwarted by the presence of dust and possibly circumstellar 
matter from nearby stellar systems. The binary progenitors 
of the helium WD channel on the other hand, although more 
abundant at most delay times, should be harder to detect as 
most of their evolutionary time is spent during the detached 
double WD phase. 

Thus far, we have only found (possibly) a weak correla- 
tion between the mass of the exploding WD and delay time, 
making it difficult to infer a connection between observed 
brightness ( 56 Ni synthesised in the explosion) and progen- 
itor age. Nevertheless, If a connection between the age of 
the primary CO WD and the production of 56 Ni can be 
made in sub-Mch explosions such that dimmer SNe la oc- 
cur among older populations, this would have very exciting 
consequences for our study. 

5.5 Conclusion 

Our standard model population synthesis indicates that 
there are potentially enough sub-Mch progenitors to ac- 
count for the rates of SNe la. Nevertheless, much uncer- 
tainty still remains regarding the formation and evolution of 
close binary stars: mass transfer and accretion efficiencies, 
effects of rotation and magnetic fields, impact of metallicity 
on stellar winds and subsequent stellar and binary evolu- 
tion, the common envelope phase, etc. Even given a large 
population of potential progenitors for sub-Mch explosions, 
there remain open questions about the explosion itself. Hy- 
drodynamical studies have previously shown that sub-Mch 
WDs with an overlying helium shell can undergo a double- 
detonation which looks like a SN la, though the real answer 
as to what fraction of these systems lead to SNe la explo- 
sions depends on specific details. Most critically, under ex- 
actly which conditions does helium ignition occur, and how 
does the nucleosynthesis proceed? 

The sub-Mch model is the first model which demon- 
strates a sufficient number of SNe la events to account for all, 
or at least some substantial fraction of, SNe la (Model Al), 
as well as two distinct formation channels with their own 
characteristic DTD: A prompt (<500 Myr) helium star 
channel originating from binaries with more massive sec- 
ondaries, and a more delayed (>500 Myr) double WD chan- 
nel originating from AM CVn-like progenitor binaries with 
lower mass. Whether some or all of the sub-Mch models ex- 
plored in this work really lead to thermonuclear explosions 
that look like normal (or some subclass of) SNe la is still a 
topic which requires further study. 
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